function f=fQR(Q1,R1)     %利润函数 
Q1x=size(Q1,2);%求Q1列数
I=size(Q1,1);%时间段数
ci=40;cI=200;r=100;
Qnew=Q1;
Rnew=R1;
      for n=1:Q1x
          Qnew(:,n)=r*(n-1)*Qnew(:,n);
          Rnew(:,n)=(n-1)*Rnew(:,n);
      end
          Rnew=ci*Rnew;
          Rnew(I,:)=cI/ci*Rnew(I,:);
      f=sum(sum(Qnew))-sum(sum(Rnew));     
end
 